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1/m Corrections for Orbitally Excited Heavy Mesons and the 1/2 — 3/2 Puzzle 
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We re-investigate the effects of the l/m c corrections on the spectrum of the lowest orbitally excited 
D-meson states. We argue that one should expect the 1 jm c corrections to induce a significant mixing 
between the two lowest lying 1 + states. We discuss the implications of this mixing and compute its 
effect on the semileptonic decays B —> D**£u and the strong D** decays. 


I. INTRODUCTION 

The spectroscopy of excited hadrons containing a 
heavy quark is determined to a large extend by the fact 
that the spin of the heavy quark decouples from the light 
degrees of freedom [1], To this end, the rotations of the 
heavy quark spin become a symmetry that is not present 
for light hadrons. As a consequence, all heavy hadrons 
(with a single heavy quark) fall into spin-symmetry dou¬ 
blets, the members of which are related by a rotation of 
the heavy quark spin. 

For the mesonic ground states, the spin symmetry dou¬ 
blets consist of the 0“ pseudo scalar meson and the 1” 
vector meson, such as (D,D*) and (B,B*). For the 
mesons with one additional unit of angular momentum 
of the light degrees of freedom two further spin symme¬ 
try doublets appear. These correspond to a total an¬ 
gular momentum j = 1/2 and j = 3/2 of the light de¬ 
grees of freedom, respectively. Coupling it to the heavy 
quark spin and taking into account the factor (—1/ +1 for 
the parity, we end up with a doublet with the quantum 
numbers (0 + ,l + ) for j = 1/2 and one with (l + ,2 + ) for 
j = 3/2. 

The same systematics continues for higher orbital ex¬ 
citations, i.e. for arbitrary £ [1], The doublets which 
emerge have the quantum numbers 

((£ — 1}(- 1 ) 4+1 ,^(- 1 ) <+1 ) for j = £—1/2, 

and 

,{£+l/~ 1)t+1 ) for j = 1+1/2. 

The decoupling of the heavy quark also has interest¬ 
ing consequences for the strong decays of excited heavy 
hadrons, since these are then governed by the light de¬ 
grees of freedom only. In this way one may predict the 
partial waves involved in the strong decays and thus one 
may predict a specific pattern of angular-momentum sup¬ 
pressions in the strong decays [1, 2], 
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In the heavy mass limit the members of the doublets 
are mass-degenerate, and the splittings between the mul- 
tiplets are mass independent. Nevertheless, these split¬ 
tings are of the order of Aqcd, since they are related to 
excitations of the light degrees of freedom. 

The l/m c corrections couple the heavy quark spin to 
the light degrees of freedom, leading also to a splitting 
within the doublets. However, in particular for the D 
mesons, the splitting between the multiplets and the 
(1/m induced) splittings within each of the multiplets 
are numerically of the same order. Thus a significant 
mixing of the two i^ 1 ' 1 states belonging to the two 
different multiplets must be expected. 

Previous analyses of the 1 /m c corrections also noticed 
the possibility of this mixing [2-6]. However, some of 
the analyses assumed that this effect is small, based on 
the observation that the two 1 + states should have dif¬ 
ferent widths, since the strong decays of the 1 + state 
with j = 3/2 correspond to p -wave transitions. Since the 
data indicate that the observed widths indeed follow this 
pattern, the mixing was assumed to be small. 

In the present paper we re-consider this effect and try 
to estimate the size of the corresponding mixing angle, 
which can be extracted from the data on the basis of 
a few model assumptions. Based on this, we extract a 
mixing angle that is larger than what has been discussed 
before, and consider the implications in particular for 
semileptonic B decays into orbitally excited D mesons. 

II. ORBITALLY EXCITED STATES AND THE 
EFFECT OF MIXING 

The orbitally excited states fall into two spin- 
symmetry doublets, classified by their angular momen¬ 
tum of the light degrees of freedom j. We shall use the 
notation 

(!£(?+))) with 

(jn*(2+))) with J = 3 / 2 - 

In the limit m c -A oo the members within each of the 
spin symmetry doublet are degenerate: 

M(D( 0+)) = M(D( 1+)) = m c + A 1/2 , 

M(D*( 1+)) = M(D*( 2+)) =m c + A 3/2 , 


( 2 ) 
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where A j is the binding energy of the mesons in the limit 

m c oo. 

Note that the splitting A 3 / 2 — ^- 1/1 between the two 
doublets does not scale with the heavy quark mass. How¬ 
ever, it is related to the binding of the light quark within 
the chromoelectric field of the heavy quark and hence 
it is of the order of Aqcd- In fact, for the D mesons, 
the current data yield a value of about 20 MeV for this 
splitting. 

Power corrections of order l/m c are induced by the 
kinetic and the chromomagnetic operator, leading to a 
Hamiltonian density of the form 

'Hi/m = 7 ^— c(iD±) 2 c+ ^-c{a- B)c. (3) 

Zm c 2m c 

In particular, the second term couples the heavy quark 
spin to the light degrees of freedom, which breaks the spin 
symmetry; however, the angular momentum of the light 
degrees of freedom j is still a good quantum number. 

We are going to consider the effects of T~Li/ m on the 
two spin symmetry doublets shown in eq. (1). First of 
all, the kinetic energy contribution only leads to a shift 
of the masses, which we shall absorb into the values of 
the masses 

Mj =m c + A j + (j), (4) 

where g^(j) I s the kinetic-energy parameter for the or- 
bitally excited meson. The second term in eq. (3), how¬ 
ever, leads to the effects, which we shall discuss in some 
detail. It is related to the chromomagnetic field induced 
by the light degrees of freedom at the location of the 
heavy quark spin. Clearly not much is known about this 
from first principles, so one has to make some assump¬ 
tions here to arrive at quantitative estimates. 

The chromomagnetic field at the location of the heavy 
quark is generated by the angular momenta of the light 
degrees of freedom, which is the orbital angular momen¬ 
tum L and the spin of the light quark s, the sum of which 
constitutes the angular momentum of the light degrees of 
freedom J. The complete angular momentum K is 

K = J + a = L + s + (?. (5) 

The key assumption we shall make is that the orbital 
angular momentum and the light quark spin have differ¬ 
ent gyro-chromomagnetic factors a' and j3', such that the 
chromomagnetic field seen by the heavy quark is not pro¬ 
portional to the total angular momentum J of the light 
degrees of freedom: 

B ~ ol L + (3's = aJ + /3s . (6) 

To this end, the second term in eq. (3) can be written as 

H 1 / m = J d 3 x-^-c(a-B)c = Pi(J-a) + P 2 (s-a), (7) 

where the operators Pi and P 2 act only on the radial 
wave functions of the light degrees of freedom. 


We note that the |D*(2 + )) and the |P(0 + )) states are 
eigenstates of the Hamiltonian, even once the above 1 jm c 
corrections are included 

H\D*{2 + )) = (m 3/2 + \g + l -g^j |D*(2+)>, (8) 

H\D(0+)) = (m 1/2 -\g + \g^j |P(0+)>, (9) 

where the mass values M 1 / 2 and M 3 / 2 are the masses de¬ 
fined in eq. (4). Furthermore, the constants g and g' are 
obtained from the radial wave functions, i.e. the matrix 
elements g ~ (Pi) and g' ~ (P 2 ); here we assume for sim¬ 
plicity that g and g' are identical for the two doublets. 
The relevant coefficients are obtained from the spin wave 
functions discussed in the appendix. 

For the two 1 + states, the second term of eq. (7) in¬ 
duces a mixing, since we have 

H\D(l + )) = (m 1/2 + \g- ±-g') |P(1+)) 

Ji (10) 

+ ^|D*(1+)), 

H\D*(l+)) = { M3/2 -lg-^-A |D*(l+)> 

72 (U) 

+ V 3 -9 , \D(1 + )), 

where we have - once more - assumed that the relevant 
matrix elements of Pi and P 2 are identified with the pa¬ 
rameters g and g' . Again, the relevant coefficients in front 
of g and g' follow from the spin wave functions given in 
the appendix. 

Clearly this analysis is drastically simplified, but an 
improvement would need the reference to some model 
for the binding dynamics of the orbitally excited mesons. 
Nevertheless, from this simplified analysis we can infer 
some information on the mixing angle from the data. 
First of all looking at the data summarized in table I 
and schematically shown in figure 1 - we note that the 
broader 1 + state has a larger mass than the narrower 1 + 
state. From heavy-quark symmetries we infer that the 
broader states decay through an S wave transition which 
is possible only for the j = 1/2 states. The j = 3/2 
states can decay (strongly) only through a D wave tran¬ 
sition which is suppressed by angular momentum, ren¬ 
dering these states narrow [1]. For this reason we make 
the assignments shown in table I. 

Starting from the heavy quark limit, the splitting 
within the two doublets is induced by l/m c effects. Ig¬ 
noring their mixing for the moment, the two 1 + states 
will cross when we switch on the 1 /to c terms, leading to 
the “level inversion” observed in the spectrum depicted 
in figure l. 1 Switching on the mixing leads to the well 


1 In fact, from the current data shown in table I only the central 
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Figure 1. Schematic of the mass hierarchy for the orbitally 
excited D mesons. 


State 

Mass [MeV] 

Width [MeV] 

j 

D( 0 + ) 

2318 ± 29 

267 ± 40 

1/2 

D(l+) 

2427 ± 40 

384 ± 120 

1/2 

D*( 1+) 

2421.4 ±0.6 

27.4 ± 2.5 

3/2 

D*( 2+) 

2462.6 ±0.6 

49 ± 1.3 

3/2 


Table I. Averages of the existing mass and width measure¬ 
ments for the D** states according to the PDG [7]. Our 
choice of the j assignment follows from total decay widths. 


known “level repulsion”, which is shown in figure 2 for the 
parameter values obtained from the fit discussed below. 
For the 1 + states this implies a mixing, which is maximal 
at the crossing point |0| = 45°. In order to explain the 
assumed “level inversion”, we must have 45° < |6*| < 90°; 
i.e., mixing beyond the crossing point. 

For a quantitative analysis, we use our simple assump¬ 
tions and encode the mixing of the two 1 + states in the 
2x2 sub matrix of the Hamiltonian 


H 


M a 
a M* 


and we define 


M = 
M* = 
A = 

a = 


M i/2 + 4# - Y 2 9 ' ’ 

5 5 

m 3/2 - 4 5 - j^g', 

M* — M = M 3/2 - M 



1/2 - 




( 12 ) 

(13) 

(14) 

(15) 

(16) 


values indicate a level inversion. Within their uncertainties, the 
masses for the two 1+ states may also allow the naively-expected 
mass hierarchy. Nevertheless, the two states are at least very 
close in mass. 


Mass splitting 
A(g)=0 s* 



Figure 2. The orbitally excited states in the j = 1/2 and 
the j = 3/2 doublets are degenerate in the infinite mass limit 
rric = 00 , corresponding to g — g' = 0. For a finite mass 
m c we find a splitting of the states. We plot the masses of 
the orbitally excited states D( 0 + ) (light-blue, long-dashed), 
D( 1 + ) (green, short-dashed), D*(l + ) (red, solid) and D*( 2 + ) 
(dark blue, dotted) as a function of the parameter g. The 
parameter g' is fixed by the condition g' = —g. The verti¬ 
cal lines indicate the best 6t value q* for q and the zero of 

A (g,g' = -g). 


For vanishing g' the Hamiltonian becomes diagonal, 
but in this case the splitting in the D*{J + ) doublet is 
twice as large as the one in the D(J + ). We note that the 
splittings within each of the doublets is equal if g = —g' 7 
which we consider as a benchmark point. This point 
is interesting to consider in terms of our simple model: 
Assuming that all overlap integrals are equal, we obtain 
from eq. (6) 

B ~ gj + g's = g{J~- s) = gL, (17) 


which would mean that the spectrum is roughly driven 
by the orbital angular momentum L of the light degrees 
of freedom. As a result, the spectrum is independent of 
the orientation of the light quark spin. However, this 
could as well be an artifact of our simple assumptions. 

To this end, the eigenstates of H (i.e., the physical 1 + 
states) are thus linear combinations of the states defined 
in the heavy-mass limit 


\D l (1+))= cos0|D(l + ))+sin0|ZT(l+)), (18) 

\D h {1+)) = — sin f? |D(1 + )) + cos 6* |-D*(1 + )), (19) 


where Dl ( Dh ) is the state with the lower (higher) eigen¬ 
value. The mixing angle 9 satisfies 


tan 26 = 


—2 a 


-4 V2g' 


A 6(Af 3 / 2 — A/ 1 / 2 ) — 9<7 + 2 g' 

and the corresponding eigenvalues are given by 
M h/l = i [(Af + M*) ± \/ A 2 + 4a 2 


( 20 ) 


(21) 
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Mixing angle 6 

9 * 



Figure 3. The mixing angle 6 as a function of g' and with 
g = —g'. The vertical line indicates the best fit point. The 
horizontal line indicates the maximum-mixing angle of 9 = 
45°. 


Mass [MeV] 

Pull [cr] 

Reference 

D{ 0+) 

2297 ± 8 

-13.75 

[8] x 

2308 ± 36 

-2.75 

[9] x 

2407 ± 41 

-0.00 

[10] 

Dl( 1+) 

2423.1 ± 1.8 

±1.24 

[H] 

2420.1 ± 0.8 

-0.95 

[12] 

2426.0 ± 3.2 

±1.61 

[13] 

Dh( 1 + ) 

2427 ± 36 

-0.62 

[9] 

2477 ± 40 

±0.69 

[14] t 

D*(2 + ) 

2462.5 ±2.7 

±0.01 

[H] 

2462.2 ±0.8 

-0.34 

[12] 

2464.5 ±2.2 

±0.92 

[10] 


Note that we have always Mr > Ml, reflecting the 
“level repulsion” of the 2x2 system. The situation of 
minimal splitting Mr —Ml corresponds to A = 0, which 
yields \9\ = 45°. This means that the contribution of 
the j = 3/2 state in the eigenvector \Dl(1 + )) starts to 
become dominant, in accordance to the observed level 
sequence. This is shown in figure 3, where the mixing 
angle 9 is plotted as a function of g' , with the constraints 
g = —g'. There, the vertical line corresponds to the best- 
fit value of g' of the fit discussed below, where we find a 
mixing angle of about 60°. 

In previous analyses (e.g. in reference [4]), a value 
for 9 of the order of ten degrees has been considered. 
However, this is not in contradiction to our discussion, 
since an interchange of the two 1 + states corresponds to 
a replacement 6* —90 — Nevertheless, we advocate 
that the magnitude of 9 is larger than considered in 
previous analyses. 

The experimental results for the masses and widths of 
the orbitally excited D mesons are presently unsettled. 
The individual results for the masses, for instance, do not 
agree well with their respective PDG averages [7]. This 
is evidenced by the application of scale factors to the 
error on the average values. Nevertheless, we attempt to 
confront our simple model with the present data. 

In a first step, we assume g' = —g, thereby reducing 
the number of model parameters to three: M 1 / 2 , M 3 / 2 , 
and g. This assumption stems from the fact, that the 
mass splitting of the doublets is equal for g = —g', as 
argued above. We then fit the three parameters to the 
four PDG averages of the measurements as listed in table 
I. Note, that we impose the additional constraint A < 0, 
based on the mass hierarchy shown in figure 1. This 
fit has one degree of freedom (d.o.f.). We reject the fit, 
since we obtain x 2 = 8.48 and a p-value of 0.36%, which 
is smaller than our a-priori threshold of 3%. While the 


Table II. Masses for the various D** measurements. We use 
up to three of the most precise mass measurements, if avail¬ 
able. We also list the pull values at the best-fit point eq. 
23. Measurements with a x mark are not included in the fit, 
see text. Uncertainties of Measurements with a j mark have 
been doubled, due to lack of estimates for the systematic un¬ 
certainty. 

masses of the two 1 + eigenstates and the 2 + state are 
modelled very well, we find a large pull of slightly less 
than 3cr for the 0 + mass. The best-fit point and the 
parameter intervals at 68% confidence level (CL) read 

M {/ 2 = (2425 ± 2)MeV, M* /2 = (2451 ± 5)MeV, 
g* = (23 ± 6)MeV. 

.( 22 ) 

In a second step, we wish to find out which of the indi¬ 
vidual measurements do not agree well with our model. 
For this purpose, we use the 11 most precise individual 
measurements, with up to three measurements per D** 
state. These measurements also enter the PDG averages, 
and they are listed in table II. Given the larger number of 
measurements, we can now lift the previous assumption 
and fit all four model parameters. We again impose the 
theory constraint A < 0. We reject this fit as well, since 
for seven d.o.f. we obtain \ 2 = 177, corresponding to a 
p-value of less than 10~ 10 . However, we observe that the 
X 2 is driven by two measurements: the measurements 
of the 0+ mass as carried out by the BaBar and Belle 
collaborations, respectively. It is therefore interesting to 
repeat this second fit without the BaBar and Belle mea¬ 
surements of 0 + masses, which we do. We find for this 
new fit the best-fit point and the 68% CL intervals 

My 2 = (2433 ± 4)MeV, M* /2 = (2450 ± 5)MeV, 

g* = (26 ± 5)MeV, g'* = (-27 ± 21)MeV. 

(23) 
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ti/ 2 (1) 

2 

Pl/2 

73/2(1) 

2 

P3/2 

GI [15] 

0.22 

0.83 

0.54 

1.50 

VD [16] 

0.13 

0.57 

0.43 

1.39 

CCCN [17, 18] 

0.06 

0.73 

0.51 

1.45 

ISGW [19] 

0.34 

1.08 

0.59 

1.76 


Table III. Results of the fit of eq. (31) to four different model 
calculations for the two form factors, (numbers taken from 
[ 20 ]). 

We also compute the goodness of fit for this reduced data 
set. We find a good fit for five d.o.f. and y 2 = 6 . 86 , 
which corresponds to a p-value of 0.23. The individual 
pull values at the best-fit point, defined via 

Mi — Mi(M^ l2 , My 2 , g*, g'*) 
pull* = - (24) 

®i 

for every measurement i, are listed in table II. Based 
on the results of this fit, eq. (23), we find approximately 
A ~ —13 MeV, as well as g' ~ —g < 0. This indicates a 
large mixing, with 9 = (+ 59 ^ 44 )° at 68 % CL. 

Our findings can be summarized as follows: Neither 
the averages, nor the three most precise measurements 
of each of the D** masses are fitted well by our simple 
assumptions. Removing two of the individual measure¬ 
ments of the 0 + mass from our analysis yields a good fit. 
In addition, most of the measurements are not in good 
agreement with each other [7]. On the basis of the current 
experimental data, the situation remains inconclusive. In 
our opinion, a simultaneous determination of the masses 
and widths of all orbitally excited states needs to be un¬ 
dertaken before a definite answer to mixing of the 1 + 
states can be given. Within such an analysis, the mixing 
effects of the decay widths could be taken into account 
as well. 


III. EFFECTS ON SEMILEPTONIC DECAYS 

Exclusive semileptonic b —> c decays are governed by 
heavy quark symmetry for both the bottom and the 
charm quark. The decays of B mesons into orbitally 
excited D mesons have been investigated in the context 
of heavy quark symmetry in [5], including also the 1 jm c 
corrections. As stated above and in contrast to [5], we 
study a scenario where the mixing is the leading l/m c 
effect. 

In the infinite mass limit, the hadronic weak transition 
currents can be described by a single Isgur Wise function 
for each multiplet. Using the spin representations for the 
states involved 


B(v) = ^[MbIs 2 ^ , 

(25) 

+ 

1 —l 

]| 

II 

T - 1 

Q 

(26) 

D;(l+-v,e) = 

(27) 

- *V) » 



one defines for some current defined by a Dirac matrix T 
the two Isgur-Wise functions as [5] 

<B(tO|&rc|£>(l + ;i/,e)> 

= 2 r 1 / 2 (W)Tr [B(v)TD(l + -, v', e)] , (28) 

(B(v)\bTc\D*(l+-v',e)) 

= V3T 3/2 (vv , )v^Tr{B(v)TD;(l+-v',e)] . (29) 

In the above the convention for the two Isgur-Wise func¬ 
tions is the same as in [5] and [21]. 


The mixing induced by l/m c effects yields for the hadronic currents 

{B\bYc\D L ) = cos 9(B\bTc\D(l + )) + sm9(B\bTc\D*(l+)), 
(B\bTc\D H ) = -sm9(B\bTc\D(l + )) + cos 9{B\bTc\D*(1 +)), 


which can be expressed in terms of the form factors Ti/ 2 (uj) and t 3 / 2 (w). 


The two form factors are parametrized [20] via 


[ 22 ]; in particular we have 


U'M 





j = 1/2 or j = 3/2 (31) 


which is expected to be a reasonable approximation, since 
the kinematic region turns out to be quite limited: 1 < 
w < 1.3. 

Not much is known about the form factors t 3 (lu). How¬ 
ever, there are sum rules constraining these form factors 


Mg ^ 9ei/ 2 |T-i/ 2 (l )| 2 , (32) 

where and go are the kinetic energy and chromo¬ 
magnet ic moment parameters, and e\j 2 is the excitation 
energy of the j = 1/2 doublet above the ground state 
doublet. Numerically, the left-hand side of this relation 
is small: g* — g 2 a <C g 2 - This motivates the so-called 
BPS limit, which leads to g 2 — /Iq = 0. In this limit we 
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Decay mode 

B (%) 

B{B~ -> D(0 + )£v) x Z3(D(0 + ) -> D+n ~) 
B(B~ -> D L (l+)lv) x B(D l ( 1 + ) -> D*+n~) 
B(B~ -> D h {l + )lv) x B{D h { 1+) -> D*+n~) 
B(B~ -»■ D*{2 + )lv) x B(D*( 2 +) -»■ D*+n~) 

0.29 ±0.05 

0.29 ± 0.14 

0.13 ±0.04 

0.078 ±0.008 


Table IV. HFAG averages for the products of the production 
branching fractions of B~ —>• D**£v decays with the decay 
branching fractions D** — > D^ + n~ [23]. The uncertainties 
arise from the squared sum of the statistical and systematical 
uncertainties. 


would have ri/ 2 (l) = 0 indicating that the theoretical 
expectation is that (in the relevant kinematic region) 

T1/2 M « t 3 / 2 (w) . (33) 

In the combined infinite mass and BPS limits this leads 
to the expectation that the decays into the D{ 1 + ) are 
heavily suppressed compared to the ones into the D* (1 + ). 
Current data do not support this, which constitutes the 
1/2 — 3/2 puzzle in semileptonic B decays. 

Various models and sum rule calculations have been 
used to obtain more information on the Tj{to). In Ta¬ 
ble III we list some of the currently used models in this 
context; note that all models reflect the relation eq. (33) 
in a more or less pronounced way. 


The differential decay rates in the infinite mass limit are well known [5]. Including the mixing eq. (30) as the leading 
l/m c effect one obtains for B —»• DlIv 


dT 

duj 


G 2 F V c 2 b M 5 B 

247T 3 


1)Voj 2 -1 


I sin 2 0(w + 1) 2 [2 (w - r L )( 1 - r L w) - (1 + r\ - 2rz / w)]|r 3 / 2 (w)| 2 
+ 2cos 2 9[2r L {u 2 - 1) + (5w - 1)(1 + r\ - 2rLw)]|r 1/2 (a;)| 2 
+ 2\/2sin0cos0(w + 1) 2 [(1 + rq,) 2 — 4 riw]i?e(r 1/ ' 2 (w) t 3/2H)} > 


while we get for B —> Dh^v 


dT G 2 F V 2 b M% 3 , ^ r 

d ^ = — 2 47r ^ r ^- 1) ^ 2 - 1 
| cos 2 9{uj + 1) 2 [2(w - r H )( 1 - r H u) - (1 +r% - 2r//w)]|r 3/2 (w)| 2 

+ 2sin 2 9[2r H (w 2 - 1) + (5w - 1)(1 + r% - 2r ff a;)]|Ti/ 2 (a;)| 2 
— 2v / 2sin0cos0(w + 1) 2 [(1 + rn) 2 — 4r//w]4?e(ri/ 2 (w)r 3 / 2 (w))|. 


(34) 


(35) 


The above rates now depend on both of the form factors as well as on the mixing angle 9. By inserting the 

expression eq. (31) with the values given in table III and using the nominal fit result cos0 ~ 0.51, we calculate the 
invidual branching ratios shown in table V. 


The first six rows of table V are the values obtained 
in the infinite mass limit, where we use the experimental 
results for the lifetimes. The last three rows are obtained 
for finite m c , where we consider mixing effects among the 
1 + states as the leading l/m c effect. 

Given the large uncertainties on the experimental in¬ 
puts, we wish to emphasize that our analysis is only 
meant as a qualitative study; i.e., to answer the ques¬ 
tion: Can mixing between the two 1 + states (at least par¬ 
tially) explain the 1/2 - 3/2 puzzle? As a consequence, 
we abstain from providing uncertainty estimates on the 
quantities in table V. Clearly we find a strong impact of 
the mixing, which roughly swaps the roles of the two 1 + 
states. We find that the inclusion of 1 /m c mixing effects 
redistributes the relative weights of the invidivual decay 
channels within the decays B~ —>• D**iv. 


Nevertheless, our results can be confronted with the 
experimental data shown in Table IV. Unfortunately, 
there are no experimental results yet available on the ab¬ 
solute branching fractions, since the branching fractions 
for the subsequent strong decays D** —> D^*)+n~ have 
not yet been measured. However, assuming that these 
subsequent decays have roughly the same branching frac¬ 
tions, the measured ratio of the two decays into 1 + states 
is approximately 


B{B~ D l { l+)tv) 
B(B~ D h (\+) lv) 


(36) 


The estimates without mixing effects (table V, row 6) 
clearly deviate from the measured ratio. On the other 
hand, our estimates with mixing effects taken into ac¬ 
count (table V, row 9) are in reasonable agreement with 
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Channel 

GI 

VD 

CCCN 

ISGW 

m c —>■ oo 

B(B~ a D{0+)£v) 

4.7- 10“ 4 

1.8 ■ 10 -4 

3.7- 10" 5 

1.0 • 10 -3 

B{B~ a D(l+)£v) 

6.4 ■ 10“ 4 

2.5 ■ 10" 4 

4.9 ■ 10" 5 

1.4- 10" 3 

B{B~ a D*(l+)£v) 

4.4 ■ 10“ 3 

2.9 ■ 10 -3 

4.0 ■ 10 -3 

4.7- 10 -3 

B(B~ a D*{2+)£v) 

7.4 ■ 10“ 3 

4.9 ■ 10 -3 

6.7- 10" 3 

8.0- 10" 3 

B(B~ a D**£v) 

1.3% 

0.82% 

1.1% 

1.5% 

B(B~ a D*(l+)£u) 

6.9 

11 

80 

3.4 

B(B~ a D(l+)£v) 

m c finite 

B(B~ a D l £v) 

CO 

1 

o 

o 

CO 

2.1 ■ 10 -3 

3.0- 10" 3 

3.2 • 10" 3 

B{B~ a D h £u) 

2.3 ■ 10“ 3 

1.3 ■ 10" 3 

1.3 ■ 10" 3 

3.1 • 10" 3 

B{B~ a D l £v) 
B(B~ a D h £v) 

1.3 

1.6 

2.3 

1.0 


Table V. Predictions for the branching fractions for the channels B —> D**£u within the various models of table III, both with 
and without mixing effects among the 1 + states. Here B(B~ A D**£v) represents the sum of four branching fractions to the 
0 + , the two 1 + and the 2 + states. Note, that for m c —> oo, we use the doublet masses M±/ 2 and M 3 / 2 as inputs, while for the 
finite m c results, we use the experimentally determined masses. 


the measurements. 


IV. EFFECTS ON THE WIDTHS OF THE 
ORBITALLY EXCITED STATES 

In [1] it has been discussed that in the infinite mass 
limit one has the relations 


^l(D*(2 + ) A Dn) oc \J^ a D 

(37) 

A(D*( 2+) a D*tt) oc ^a D 

(38) 

A(D*( 1+) a D*tt) oc a D 

(39) 

A{D*(1+) a£>7t)=0 

(40) 

A(D(l + ) a D*tt) oc as 

(41) 

A{D{1+) A Dtt) = 0 

(42) 

H(I?(0 + ) A Dn) oc as 

(43) 

A{D( 0+) A D* tt ) = 0 

(44) 


where an and as are the amplitudes for the D -wave and 
the S'-wave decays of the D** states. 

We assume that these modes dominate the total 
widths, such that 

r tot (D**) = r(D** -A Dtt) + T(D** A D*ir ), (45) 

so we obtain the predictions in the heavy quark limit 

r tot (D*(2+)) = M 2 = r tot (D*(l + )), 
r tot (£>(o + )) = M 2 = r tot (D(i+)), 

where we ignore small phase space differences of the order 
of twenty percent. Since we expect the D -wave amplitude 


to be suppressed relative to the S- wave amplitude by the 
usual angular momentum factors, we arrive at the well 
known conclusion that the j = 1/2-doublet states are 
broader than the states of the j = 3/2 doublet. 

Including the mixing induced by the 1 / m c terms yields 
the relations 

A{D L {1+) a Dtt) = 0 = A{D h { 1+) A Dtt) , (47) 

A(Dl(1 + ) a D*tt ) = as cos 9 + an sin 0 , (48) 

A(Dh{1 + ) A D*n) = an cos 9 — as sin# . (49) 

Within the differential decay width, interference terms 
between the S and D wave amplitudes arise. However, 
they drop out after integrating over the D*tt helicity an¬ 
gle. Consequently we have 

T t o t (D L ) = r L 

~ T(Dl —> D*tt) = |a5| 2 cos 2 9 + |a£>| 2 sin 2 6 , (50) 

and 

T tot (D H ) = T h 

~ T(Dh —> D*tt) = \au\ 2 cos 2 9 + |as| 2 sin 2 9 , (51) 

where we again ignore small phase space differences. 

The experimental situation on strong decays of the D** 
states is shown in table I and is not yet conclusive. Nev¬ 
ertheless, we can get a qualitative picture by assuming 
that we can extract |as| and \au\ from the widths of the 
0 + and the 2 + states, respectively, and insert these into 
eq. (50). From this we obtain 

T l ~ 240 MeV, T H ~ 80 MeV, 

where we again do not consider the experimental uncer¬ 
tainties, since we only aim at the qualitative picture. 




Thus the observed pattern is not in contradiction to 
a large mixing of the 1 + states, although the observed 
factor of about two between the widths of the two narrow 
states still remains unexplained. 


V. CONCLUSION 

We have discussed a scenario where the mixing of 
the two 1 + states of the first orbitally excited charmed 
mesons is assumed to be large. Such a large mixing can 
still be accommodated with the data. It is supported by 
very simple arguments on the physics origin of such a 
mixing and some assumptions on the size of some ma¬ 
trix elements. However, the input for any estimate is the 
data on the masses and the widths. The data are not 
yet conclusive, at least not for the broader one of the 1 + 
states. 

If such a large mixing is indeed present, it will also 
have consequences for the view on spectroscopy from the 
perspective of the heavy quark limit. At least for the 
charm mesons this means that the spin-symmetry dou¬ 
blets will have a significant mixing for the states with the 
same quantum numbers, which is induced by interactions 
that are formally of subleading order but are numerically 
significant. 

For the semileptonic decays of B mesons such a mixing 
would soften the 1/2 - 3/2 puzzle at least for the the In¬ 
states, since a mixing with a significant angle will reduce 
the difference of the two rates. Nevertheless, in order 
to pin down if the mixing is really the solution to this 
puzzle, more data on the semileptonic decays as well as 
on the masses and the widths of the orbitally excited D 
meson states will be required. 
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Appendix A: Spin Wave Functions 

In this appendix we give the explicit formulae for the 
coupling of the angular momentum L = 1 and the light 
quark spin s = 1/2. The spin wave functions for the case 
j = 1/2 read 

\j ~ 1/2, +1/2) = y||l)|-l/2), - y||0)|l/2)i (Al) 

\j = 1/2, 1/2) = sj\ |0)|-l/2), - y||-l)|l/2), (A2) 
and for j = 3/2 we get 

\j = 3/2,+3/2) = |+ 1 )|+ 1 / 2 ); (A3) 

I j = 3/2,+ 1 / 2 ) = ^1)1-1/2), + y||0)|l/2), (A4) 

I j = 3/2,- 1 / 2 ) = y||0)|-l/2), + ^1-1)11/2); (A5) 
|i = 3/2, -3/2) = | — 1 ) | — l/ 2 ) ; (A 6 ) 

where the first ket vector is for the angular momentum, 
while \-}i denotes the spin of the light quark. 

The above states have to be combined with the heavy 
quark spin in order to obtain the spin wave functions of 
the D** mesons. Since we are interested in the mixing of 
the 1 + states, we concentrate on these states and obtain 


|D(1+), M = 1) = | j = 1/2,1/2) 11/2) jy (A7) 

\D(1 + ), M = 0) = ^ (V = 1/2,1/2)| - 1/2) h + |j = 1/2, —1/2)11/2)^ (A8) 

\D(1+),M = -1) = | j = 1/2,-1/2)| - 1/2) H (A9) 

and 

\D*(1+),M = 1) = y|| j = 3/2,3/2)| - 1/2) H ~\j\\j = 3/2,1/2)11/2)^ (A10) 

\D*(1+),M = 0) = ^1 j = 3/2,1/2)| - 1/2) H -\j\\3= 3/2, - l/2)|l/2) ff (All) 

|D*(1+), M = -1) = ^1 j = 3/2, -1/2)| - 1/2) h - ^1 j = 3/2, -3/2)|l/2) H (A12) 


Combining this with eqs. (A1-A6) yields the wave func- tions which can be used to discuss the mixing. 
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The spin-spin coupling can be computed using 


which yields for any M the result 


(s- a) = -(s+er_ + s_<7 + ) + s 3 cr 3 


(s • a)\D(l+)) = -1 1£>(1+)) + ^V(l + )> , (A14) 
(s • a)\D* (1 + )) = -^|D*(l + )> + ^ \D(1+ )). (A15) 

(A13) 
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